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ABSTRACT 

It is shown that the finite speed of gravity affects very-long baseline interferometric 
observations of quasars during the time of their line-of-sight close angular encounter with 
Jupiter. The next such event will take place in 2002, September 8. The present Letter 
suggests a new experimental test of general relativity in which the effect of propagation 
of gravity can be directly measured by very-long baseline interferometry as an excess 
time delay in addition to the logarithmic Shapiro time delay (Shapiro, I. I., 1964, Phys. 
Rev. Lett., 13, 789). 

Subject headings: gravitation - relativity - techniques: interferometric - (galaxies:) 
quasars: individual (QSO J0842+1835) 

1. Introduction 

Experimental verifications of basic principles underlying Einstein's general relativity theory 
are important for fundamental physics. Numerous tests of general relativity in the solar system 
(Will 1993) and in the binary pulsar PSR 1913+16 (Taylor 1994) confirm its validity up to the 
precision of 1% or slightly better. 

It is worth emphasizing that all of the Solar system tests of general relativity have relied 
upon the Schwarzschild solution and could say nothing about effects of retardation associated with 
the propagation of gravity. Einstein's theory of general relativity predicts that if the second time 
derivative of the quadrupole moment of a gravitating system (e.g., a binary pulsar) is not zero, 
the system emits gravitational waves that travel outward at the speed of light. Indirect evidence 
for the existence of gravitational waves consists of observations of the inspiraling orbits of binary 
pulsars (Taylor 1994). Observations of binary pulsar inspiral for PSR 1913+16 agree at the level of 
~ 1% with predictions based on the emission of energy by the binary in the form of (quadrupole) 
gravitational radiation (Kopeikin 1985; Schafer 1985; Damour 1987). This also provides an indirect 
evidence that gravity waves must travel at the speed of light (Damour 1987; Will 1993). However, 
no one has yet directly detected gravitational waves, let alone measured their speed. 
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It is the purpose of this Letter to point out that observing propagation of light through the 
gravitational field of the solar system can serve as a tool for measuring effects associated with the 
finite speed of propagation of gravity. This is based on our previous papers (Kopeikin 1997; Kopeikin 
et al. 1999; Kopeikin Sz Schafer 1999), where we have developed a post-Minkowskian approach for 
solving the problem of propagation of light through time-dependent gravitational fields in the 
geometric optics approximation. This approach is based on making use of the retarded Lienard- 
Wiechert-type solutions of the linearized Einstein equations and allows one to find a smooth analytic 
representation of light-ray trajectory for arbitrary locations of the source of light and observer 
without imposing any restrictions on the motion of the light-ray deflecting bodies. We find that 
electromagnetic signals interact with gravitating bodies only through the retarded gravitational 
field of the bodies. That is, if a light-ray deflecting body moves with respect to a chosen coordinate 
system the temporal variation in its gravity field must take time to reach the electromagnetic signal 
in order to perturb its trajectory. This observation constitutes the main idea of the proposed VLBI 
test of the propagation of gravity elaborated in the following sections of the present Letter. 

In 1964, I. I. Shapiro (1964) suggested that the gravitational deflection of light by the Sun - 
one of the three classical effects of general relativity analyzed by Einstein — could be measured more 
accurately at radio wavelengths with interferometric techniques than at visible wavelengths with 
the available optical techniques. His idea led to stringent experimental limitations on the parameter 
7 of the parametrized post-Newtonian formalism (Shapiro 1967), thereby strongly restricting the 
number of viable theories of gravity. Our post-Minkowskian approach gives in the case of static 
spherically-symmetric field the same result as predicted by Shapiro. Furthermore, we are able to 
calculate additional corrections to the Shapiro time delay related to the non-stationarity of the 
gravitational field of the solar system and associate them with the finite speed of propagation of 
gravity. The largest contribution to the non-stationarity of gravitational field of the solar system 
comes about via the orbital motions of the most massive planets — Jupiter and Saturn. Therefore, 
it is reasonable to undertake an attempt to detect the effect of propagation of gravity by observing 
very accurately the deflection of light rays from a background source of light (quasar) caused by 
the motion of Jupiter or Saturn. This is the essence of our proposed VLBI test of general relativity 
discussed in the present paper. 

In section 2 we consider the basic formula for relativistic time delay in time-dependent gravi- 
tational fields. Section 3 outlines the basic principles of measurement of the effect of propagation 
of gravity in radio interferometric experiments. Section 4 is dedicated to the description of the 
proposed VLBI experiment and gives numerical estimates for the Shapiro time delay in the field of 
Jupiter and for the effect of propagation of gravity. Finally, in section 5 we discuss the proposed 
VLBI experiment. 
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2. Relativistic time delay in time-dependent gravitational fields 

Let us assume that the gravitational field is generated only by the Solar system's bodies and 
there is a global four-dimensional coordinate system (t, x) with origin at the barycenter of the Solar 
system. The total time of propagation of an electromagnetic signal from the point (to, Xq) (quasar) 
to the point (t, x) (observer) is given by the expression 



t = -\x 

c 



x \ + A(t,*o) • 



(1) 



Here, \x — x \ is the usual Euclidean distance between the points of emission, x , and observation, 
x, of the photon, and A(t, to) is the relativistic time delay produced by the gravitational field of 
the moving gravitating bodies of the Solar system. The function A(i, to) is given by (Kopeikin & 
Schafer 1999) 
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where m a is the mass of the ath body, t* is the time of the closest approach of electromagnetic signal 
to the barycenter of the Solar system 1 , x a (t) are coordinates of the ath body, v a (t) = dx a (t)/dt is 
the (non-constant) velocity of the ath light-ray deflecting body, k is the unit vector from the point 
of emission to the point of observation, s is a retarded time obtained by solving the gravitational 
null cone equation for the time of observation of photon t, and so is found by solving the same 
equation written down for the time of emission of the photon t 



s + -\x - x a (s)\ = t , 
c 



s + - \x - £c a (s )| = t . 
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(3) 



It is remarkable that formula (2) does not impose any restriction on the motion of bodies since 
we did not use any explicit solutions of their equations of motion while calculating the light ray 
propagation. The formula (2) can be further simplified by performing integration by parts which 
yields 2 
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where r a = x - x a (s), r 0a = x - x a (s ), r a = \r a \, r 0a = \r 0a \, v a = x a (s), v a0 = x a (s ), a a = 
v a (s), and the retarded times s and so are calculated from the gravitational null-cone equations 
(3). 



lr The time t* is used in calculations as a mathematical tool only. 
2 After integration we omitted for simplicity terms of order Va/c 2 . 
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Our expression (4) vastly extends the region of applicability of Shapiro's work (Shapiro 1967) 
for it is valid for the case of arbitrary-moving bodies whereas the calculations by Shapiro had been 
done only for the case of static Schwarzschild field. Earlier work on improving the Schwarzschild 
approximation in the case of light passing through the gravitational field of moving bodies should 
be noted (Hellings 1986; Klioner 1991; Klioner Sz Kopeikin 1992; Sovers, Fanselow & Jacobs 1998); 
however, our expression (4) is much more general and completely solves the problem. We empha- 
size that solution (4) is valid everywhere both inside and outside of the Solar system at arbitrary 
distances including its near, intermediate, and wave zones, and, for this reason, represents a smooth 
analytic solution of the equations of light-ray geodesies for arbitrarily located source of light, ob- 
server, and the system of the light-ray deflecting bodies. However, the most important observation 
is that our formula (4) allows one to single out the retardation effects uniquely associated with the 
finite speed of propagation of gravity and indicates that for the correct calculation of the relativis- 
tic time delay positions of the gravitating bodies must be taken at the retarded times so an d s 
corresponding to the instants of emission to an d observation t of the electromagnetic signal 



3. Principles of Measurement of the Propagation of Gravity by VLBI 

Very Long Baseline Interferometry (VLBI) has recently reached a precision in measurements 
of differential phase delay of order 10 -12 second (ps), which makes it presently the most accurate 
technique for measuring relativistic time delay (Beasley & Conway 1995; Honma, Kawaguchi & 
Sasao 2000). In order to discuss how to measure the relativistic effect of propagation gravity by 
VLBI one transforms the formula (4) to the form used in VLBI analysis and called the differential 
VLBI time delay (Walter & Sovers 2000). In doing this conversion, we note that formula (4) is 
also valid in scalar-tensor theories after replacement of the universal gravitational constant G — > 
G(l + 7)/2, where 7 is one of the parameters of the PPN formalism (Will 1993). 

Let us assume now that there are two earth-based VLBI stations and the plane front of elec- 
tromagnetic waves from a quasar propagates towards the Earth. Taking two rays from the wave 
front and subtracting equation (1) for the first light ray from that for the second ray yields 

1 1 

*2 - h = -\x2 - x \ - -\xi - x \ + A(£i,t 2 ) , (5) 

where A(ii, i 2 ) = A(<2, to) — A(ti, to). Here t\ and i 2 are coordinate times of arrival of the electro- 
magnetic signal from the source of electromagnetic waves (quasar) to the first and second stations, 
respectively, and X\ and X2 are spatial coordinates of the first and second VLBI stations with 
respect to the barycentric frame of the solar system that is chosen as the primary non-rotating 
reference frame (Walter & Sovers 2000). 

The difference A(ii,t 2 ) is obtained from the main formula (4) after long and tedious calcula- 
tions and analysis of residual terms. It can be proved (Kopeikin & Schafer 1999) that with accuracy 
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better than 1 ps, 

Afe.t.) - (l +T )£;g£(l + g*(S)) ln r 'f'i+^'f'; , (6) 

V c / ^2a(s 2 ) + K-r 2a (s 2 ) 

where is the unit vector from the barycenter of the solar system to the quasar; v a (s\) is the 
velocity of the ath gravitating body at time si, r± a = \ri a \, r 2a = |r 2a |, ^ia(si) = xi(ti) — x a (si) 
and r 2a (s 2 ) = x 2 (t 2 ) — x a{ s 2) \ moreover, the retarded times s\ and s 2 are calculated according to 

si = h \xi(ti) - x n (si)| , s 2 = *2 \x 2 {t 2 ) - x a (s 2 )\ ■ (7) 

c c 

The effect of propagation of gravity appears in equation (6) as a displacement of the light-ray 
deflecting bodies from their present to retarded positions. We note that velocities of the gravitating 
bodies v a are small with respect to the speed of gravity and the time taken by the light ray to reach 
the observer after passing by the gravitating body is much smaller than its orbital period. Thus, 
one can expand the retarded positions x a (si) (i = 1,2) of the bodies in Taylor series around their 
present positions x a (tj) taken at the arrival times £j (i = 1, 2). Neglecting all terms quadratic with 
respect to velocities of the gravitating bodies and/or proportional to their accelerations one gets 
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where B = x 2 (t\) — X\(t\) is the baseline between the two VLBI stations, Vi a = Xj(tj) — x a (ti) is 
the difference of coordinates of the ith VLBI station and ath gravitating body taken at the time of 
arrival of the radio signal to the ith station, iV la = ri a /Vi , v a = v a (ti) and all quantities in the 
second term on the right side of equation (8) are also evaluated at time t±. 

The logarithmic term on the right side of equation (8) was first discovered by Shapiro (1967). 
However, the finite speed of propagation of gravity leads to the appearance of an additional (second) 
term on the right side of equation (8) as well. Had the speed of gravity been set equal to infinity 
as it is done in Newtonian gravitation, the second term on the right side of equation (8) would be 
identically equal to zero. Experimental confirmation of the existence of such a term would directly 
prove that gravity propagates with a finite speed. 



4. The Proposed Radio Interferometric Experiment 

In 2002, on September 8 Jupiter will pass at the angular distance of about 3.7 arcminutes 
from the quasar QSO J0842+1835, making an ideal celestial configuration for measuring the speed 
of propagation of gravity by using the differential VLBI technique. During the passage of Jupiter 
near the quasar line of sight, the time-dependent impact parameter |£| of the light ray from the 
quasar with respect to Jupiter will be always small as compared with the Earth-Jupiter distance, 
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r&j, which is approximately equal to 6 astronomical units. Hence, it is convenient to introduce 
the unit vector n along the direction of the impact parameter according to the definition £ = |£| n 
such that 

Nu = -cos 6 K + sin9 n , (9) 

where the subscript J refers to Jupiter and 6 ~ |£|/r©j is the (small) angle between the undisturbed 
geometric positions of the quasar and Jupiter. The baseline B during the time of the experiment 
will be much smaller than |£|; therefore, the following expansions are valid 

N u = -(l- d -pjK + 9n + 0(9 3 ) , (10) 
r 2a + K-r 2a = r la + K-r la + N la -B + K-B + 0(B 2 ) . (11) 



Other than Jupiter, the Earth and the Sun also contribute significantly to the gravitational 
time delay and must be included in the data processing algorithm in order to extract the effect 
of propagation of gravity unambiguously and to measure its speed. Accounting for this fact and 
making use of the expansions (10)-(11), we can recast formula (8) into the following convenient 
form 



A(ti,* 2 ) = A e + A -(l + 7) 



2GMj n-B 
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(12) 



where we have introduced a new phenomenological parameter 5 parameterizing the effect of propa- 
gation of gravity in data processing algorithms. Thus, we emphasize that there are two relativistic 
parameters to be measured in the VLBI experiment in order to test the validity of general relativity 
theory — the PPN parameter 7 and the gravity propagation parameter 5. The best experimental 
measurement of the parameter 7 is due to Lebach et al. (1995), who obtained 7 = 0.9996 ± 0.0017 
in excellent agreement with general relativity. The primary goal of the new experimental test of 
general relativity proposed in the present Letter is to set direct observational limits on the parame- 
ter S which will measure the effect of retardation in propagation of gravity by the moving Jupiter. 
According to the Einstein theory of relativity one must expect that the numerical value of the 
parameter 5 must be equal to zero. 

We emphasize that our consideration is fully based on the Einstein theory of general relativ- 
ity in which the speed of propagation of gravity is equal to the speed of light in vacuum. Will 
(1971) analyzed the case when the two speeds are not equal and showed that the resulting theory 
would have a non-zero value of the PPN parameter a 2 = (c/c g ) 2 — 1, where c g is the value of the 
propagation speed of gravity in the rest frame of the universe. 3 This parameter has been strongly 
bounded by consideration of various "preferred frame effects" . The present best limit (Nordtvedt 
1987) is \a 2 \ < 4 x 10~ 7 . The experiment proposed in the present Letter is not designed to compete 



3 We would like to emphasize that our parameter 5 is just a fiting parameter and has no any relation to the PPN 
parameter Q2 as we work in the framework of general relativity where c g = c. 
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with this tight bound, rather its main purpose is to measure the retardation effect (parameter S) 
associated with finite speed of propagation of gravity which was never observed before. 

The first term on the right side of equation (12) describes the Shapiro time delay due to the 
gravitational field of the Earth. It can reach 21 ps for a baseline of B ~ 6000 km. The second term 
on the right hand side of equation (12) describes the Shapiro time delay due to the Sun. It can 
vary (for B ~ 6000 km) from 17 x 10 4 ps for the light ray grazing the Sun's limb to only 17 ps when 
the direction to the source of light is opposite to the Sun (Klioner 1991). The third term on the 
right side of equation (12) is the standard Shapiro time delay due to the static gravitational field of 
Jupiter. Finally, the fourth term on the right side of equation (12) is the time delay caused by the 
effect of propagation of gravity. Numerical estimates of the Shapiro time delay Aj and the gravity 
propagation time delay Ajp caused by Jupiter can be obtained directly from equation (12). One 
finds for 7 = 1 and 5 = 0, 

where a© = 1.75" is the relativistic deflection of light by the Sun for the case of light grazing the 
Sun's limb, R & is the radius of the Sun, Rj is the radius of Jupiter, i?® is the radius of the Earth, 
9 j = 0.27' is the visible angular radius of Jupiter at the distance of r©j = 6 AU, 6 = 3.7' is the 
minimum angular distance from Jupiter to the quasar in 2002, on September 8, and vj is the orbital 
speed of Jupiter. 

It is worth emphasizing that Aj ~ while the gravity propagation time delay Ajp ~ 9~ 2 . 
For this reason, the effect of propagation of gravity is not strongly suppressed by the presence 
in Ajp of the small factor vj/c = 4.5 x 10~ 5 , and, it turns out to be measurable using VLBI 
techniques. Indeed, using the numerical values for the parameters in equation (13) for the case of 
the Jupiter-quasar "encounter" in 2002, September 8, we obtain (B ~ i? ffi ) 

aj ~ 1.26 mas , Aj ~ 137 ps , Ajp ~ 6 ps , (14) 

so that Ajp/Aj ~ 0.04. Hence, the correction due to the propagation of gravity would be ~ 4% of 
the Shapiro time delay. With the present differential VLBI accuracy it appears that the new effect 
(Ajp ~ 6 ps) is measurable with careful observations. Indeed, one can achieve a phase-referencing 
VLBI accuracy of a few degrees of phase (averaged over a 10- hour observation) and can detect 
relative astrometric position changes of less than 0.010 mas (~ 1 ps) with the existing VLBA and 
Effelsberg radio antennas (E. Fomalont, private communication). 

Formula (13) also explains why it is not so effective to observe the effect of the propagation 
of gravity in the field of the Sun. It turns out that for the case of the Sun the angle 6 can not be 
made as small as in the case of Jupiter, because the Sun is closer to the Earth and, in addition, its 
radius is ten times larger than that of Jupiter. This makes the effect of the propagation of gravity 
in the field of the Sun < 10 ps for the light ray grazing the Sun's limb. 
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5. Discussion of the Proposed VLBI Experiment 



The effect of propagation of gravity appears in addition to the logarithmic Shapiro time delay as 
an excess delay of ~ 4%. The Shapiro time delay Aj caused by the static part of the gravitational 
field of Jupiter was first measured in 1991 (Treuhaft and Lowe 1991). Klioner (1991) used the 
post-Newtonian expression for the metric tensor under the assumption that the planets of the solar 
system move uniformly along straight lines. He obtained the same result for Ajp formally; however, 
due to the specific assumptions used for the derivation of the near field post-Newtonian metric, 
he was not able to interpret it properly as the effect caused by the finite speed of propagation 
of gravity. This is because there was no proof that the solution of the light geodesic equations 
obtained on the basis of the near field post-Newtonian metric can be smoothly matched with the 
far field metric describing emission gravitational waves by the orbital motion of planets in the Solar 
system. The present Letter provides evidence that this near field post-Newtonian metric directly 
depends on the speed of propagation of gravity in general relativity. 

Other authors (Hellings 1986; Klioner 1991; Sovers, Fanselow & Jacobs 1998) have also at- 
tempted to derive more exact expressions for the Shapiro effect in the case of time-dependent 
gravitational field of the Solar system. These authors used the post-Newtonian metric for the cal- 
culation of the Shapiro effect with positions of bodies taken at the time of the closest approach t* 
of a radio signal to the body 



The instant of time t* a is numerically close to the retarded time s calculated as the solution of 
the gravitational null cone equation (3) and, for this reason, makes a good approximation for the 
calculation of the Shapiro time delay. However, in the general case when the impact parameter is 
large, t* and s can have very different numerical values. Indeed, the physical meanings of the time 
of the closest approach t* and the retarded time s are crucially different. As it is seen from its 
definition (15), the time t* is calculated only approximately using the finite speed of propagation of 
light, while the retarded time s emerges in our more exact post-Minkowskian analytic calculations 
(4) due to the finite speed of propagation of gravity. It is because of this significant difference that 
the VLBI experiment in 2002, September 8 will probe the non-trivial effect of the propagation of 
gravity. 

It is also useful to see how the effect of the propagation of gravity can be explained qualitatively 
without doing lengthy calculations. For this purpose, we note that equation (6) indicates that the 
Shapiro time delay is actually a function of the retarded time Aj(s). Expanding this function in 
Taylor series around the time of observation t, one finds 



t *a = t ~ ~ k - ( x ~ x a) ■ 



(15) 



Aj(s) = Aj(t) + Aj{t)(s - t) + ... 



a Mt)+Mt) (|)(iK) 



+ ... = Aj(t) + A JP (t) + ... , 



(16) 



where Aj(t) ~ —(9/9)Aj using equation (13), s — 



t ~ — 



r^j/c using equation (3), and Ajp is given 
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by the approximate expression (13). 

One can see that the delay Ajp due to the finite speed of propagation of gravity can become 
comparable with the standard (instantaneous) prediction for the Shapiro time delay Aj(i) in case 
9 ~ vj/c. This can never happen for Jupiter and the other solar system bodies, but can be 
interesting for VLBI experiments if one were able to find a binary star close to the line of sight 
of a quasar. Orbital motion of the binary star would cause periodic modulation of the Shapiro 
delay in the time of propagation of light from the quasar to the observer, which can be used for the 
determination of the effect of propagation of gravity outside the solar system (Kopeikin et al. 1999; 
Kopeikin &; Gwinn 2000). It is worth noting that the effect of propagation of gravity can also be 
studied in the timing of binary pulsars with nearly edgewise orbits and by the doppler tracking of 
spacecraft in deep space. 

We are indebted to Prof. B. Mashhoon who participated in numerous scientific discussions, 
carefully read the manuscript and made many essential suggestions for its improvement. We are 
grateful to Dr. S.M. Kudryavtsev for calculating the date of the Jupiter's encounter with QSO 
J0842+1835. Thanks are also due to Dr. E. Fomalont, Dr. L. Petrov, Prof. K. Nordtvedt, Prof. 
J. Burns, Prof. C. R. Gwinn, and Prof. O. Sovers for interesting and stimulating discussions and 
suggestions. 
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